Long-read sequencing and optical mapping generates near T2T assemblies that resolves a centromeric translocation.

Long-read genome sequencing (lrGS) is a promising method in genetic diagnostics. Here we investigate the potential of lrGS to detect a disease-associated chromosomal translocation between 17p13 and the 19 centromere. We constructed two sets of phased and non-phased de novo assemblies; (i) based on lrGS only and (ii) hybrid assemblies combining lrGS with optical mapping using lrGS reads with a median coverage of 34X. Variant calling detected both structural variants (SVs) and small variants and the accuracy of the small variant calling was compared with those called with short-read genome sequencing (srGS). The de novo and hybrid assemblies had high quality and contiguity with N50 of 62.85 Mb, enabling a near telomere to telomere assembly with less than a 100 contigs per haplotype. Notably, we successfully identified the centromeric breakpoint of the translocation. A concordance of 92% was observed when comparing small variant calling between srGS and lrGS. In summary, our findings underscore the remarkable potential of lrGS as a comprehensive and accurate solution for the analysis of SVs and small variants. Thus, lrGS could replace a large battery of genetic tests that were used for the diagnosis of a single symptomatic translocation carrier, highlighting the potential of lrGS in the realm of digital karyotyping.

Exonic trinucleotide repeat expansions in ZFHX3 cause spinocerebellar ataxia type 4: A poly-glycine disease.

Autosomal-dominant ataxia with sensory and autonomic neuropathy is a highly specific combined phenotype that we described in two Swedish kindreds in 2014; its genetic cause had remained unknown. Here, we report the discovery of exonic GGC trinucleotide repeat expansions, encoding poly-glycine, in zinc finger homeobox 3 (ZFHX3) in these families. The expansions were identified in whole-genome datasets within genomic segments that all affected family members shared. Non-expanded alleles carried one or more interruptions within the repeat. We also found ZFHX3 repeat expansions in three additional families, all from the region of Skåne in southern Sweden. Individuals with expanded repeats developed balance and gait disturbances at 15 to 60 years of age and had sensory neuropathy and slow saccades. Anticipation was observed in all families and correlated with different repeat lengths determined through long-read sequencing in two family members. The most severely affected individuals had marked autonomic dysfunction, with severe orthostatism as the most disabling clinical feature. Neuropathology revealed p62-positive intracytoplasmic and intranuclear inclusions in neurons of the central and enteric nervous system, as well as alpha-synuclein positivity. ZFHX3 is located within the 16q22 locus, to which spinocerebellar ataxia type 4 (SCA4) repeatedly had been mapped; the clinical phenotype in our families corresponded well with the unique phenotype described in SCA4, and the original SCA4 kindred originated from Sweden. ZFHX3 has known functions in neuronal development and differentiation n both the central and peripheral nervous system. Our findings demonstrate that SCA4 is caused by repeat expansions in ZFHX3.

Complete Mitochondrial DNA Genome Variation in the Swedish Population.

The development of complete mitochondrial genome (mitogenome) reference data for inclusion in publicly available population databases is currently underway, and the generation of more high-quality mitogenomes will only enhance the statistical power of this forensically useful locus. To characterize mitogenome variation in Sweden, the mitochondrial DNA (mtDNA) reads from the SweGen whole genome sequencing (WGS) dataset were analyzed. To overcome the interference from low-frequency nuclear mtDNA segments (NUMTs), a 10% variant frequency threshold was applied for the analysis. In total, 934 forensic-quality mitogenome haplotypes were characterized. Almost 45% of the SweGen haplotypes belonged to haplogroup H. Nearly all mitogenome haplotypes (99.1%) were assigned to European haplogroups, which was expected based on previous mtDNA studies of the Swedish population. There were signature northern Swedish and Finnish haplogroups observed in the dataset (e.g., U5b1, W1a), consistent with the nuclear DNA analyses of the SweGen data. The complete mitogenome analysis resulted in high haplotype diversity (0.9996) with a random match probability of 0.15%. Overall, the SweGen mitogenomes provide a large mtDNA reference dataset for the Swedish population and also contribute to the effort to estimate global mitogenome haplotype frequencies.

The long non-coding RNA LINC00707 interacts with Smad proteins to regulate TGFß signaling and cancer cell invasion.

BACKGROUND: Long non-coding RNAs (lncRNAs) regulate cellular processes by interacting with RNAs or proteins. Transforming growth factor ß (TGFß) signaling via Smad proteins regulates gene networks that control diverse biological processes, including cancer cell migration. LncRNAs have emerged as TGFß targets, yet, their mechanism of action and biological role in cancer remain poorly understood. METHODS: Whole-genome transcriptomics identified lncRNA genes regulated by TGFß. Protein kinase inhibitors and RNA-silencing, in combination with cDNA cloning, provided loss- and gain-of-function analyses. Cancer cell-based assays coupled to RNA-immunoprecipitation, chromatin isolation by RNA purification and protein screening sought mechanistic evidence. Functional validation of TGFß-regulated lncRNAs was based on new transcriptomics and by combining RNAscope with immunohistochemical analysis in tumor tissue. RESULTS: Transcriptomics of TGFß signaling responses revealed down-regulation of the predominantly cytoplasmic long intergenic non-protein coding RNA 707 (LINC00707). Expression of LINC00707 required Smad and mitogen-activated protein kinase inputs. By limiting the binding of Krüppel-like factor 6 to the LINC00707 promoter, TGFß led to LINC00707 repression. Functionally, LINC00707 suppressed cancer cell invasion, as well as key fibrogenic and pro-mesenchymal responses to TGFß, as also attested by RNA-sequencing analysis. LINC00707 also suppressed Smad-dependent signaling. Mechanistically, LINC00707 interacted with and retained Smad proteins in the cytoplasm. Upon TGFß stimulation, LINC00707 dissociated from the Smad complex, which allowed Smad accumulation in the nucleus. In vivo, LINC00707 expression was negatively correlated with Smad2 activation in tumor tissues. CONCLUSIONS: LINC00707 interacts with Smad proteins and limits the output of TGFß signaling, which decreases LINC00707 expression, thus favoring cancer cell invasion. Video Abstract.

Copy number variations and their effect on the plasma proteome.

Structural variations, including copy number variations (CNVs), affect around 20 million bases in the human genome and are common causes of rare conditions. CNVs are rarely investigated in complex disease research because most CNVs are not targeted on the genotyping arrays or the reference panels for genetic imputation. In this study, we characterize CNVs in a Swedish cohort (N = 1,021) using short-read whole-genome sequencing (WGS) and use long-read WGS for validation in a subcohort (N = 15), and explore their effect on 438 plasma proteins. We detected 184,182 polymorphic CNVs and identified 15 CNVs to be associated with 16 proteins (P < 8.22×10-10). Of these, 5 CNVs could be perfectly validated using long-read sequencing, including a CNV which was associated with measurements of the osteoclast-associated immunoglobulin-like receptor (OSCAR) and located upstream of OSCAR, a gene important for bone health. Two other CNVs were identified to be clusters of many short repetitive elements and another represented a complex rearrangement including an inversion. Our findings provide insights into the structure of common CNVs and their effects on the plasma proteome, and highlights the importance of investigating common CNVs, also in relation to complex diseases.

Genomic, transcriptomic and epigenomic sequencing data of the B-cell leukemia cell line REH.

OBJECTIVES: The aim of this data paper is to describe a collection of 33 genomic, transcriptomic and epigenomic sequencing datasets of the B-cell acute lymphoblastic leukemia (ALL) cell line REH. REH is one of the most frequently used cell lines for functional studies of pediatric ALL, and these data provide a multi-faceted characterization of its molecular features. The datasets described herein, generated with short- and long-read sequencing technologies, can both provide insights into the complex aberrant karyotype of REH, and be used as reference datasets for sequencing data quality assessment or for methods development. DATA DESCRIPTION: This paper describes 33 datasets corresponding to 867 gigabases of raw sequencing data generated from the REH cell line. These datasets include five different approaches for whole genome sequencing (WGS) on four sequencing platforms, two RNA sequencing (RNA-seq) techniques on two different sequencing platforms, DNA methylation sequencing, and single-cell ATAC-sequencing.

Long-read whole-genome analysis of human single cells.

Long-read sequencing has dramatically increased our understanding of human genome variation. Here, we demonstrate that long-read technology can give new insights into the genomic architecture of individual cells. Clonally expanded CD8+ T-cells from a human donor were subjected to droplet-based multiple displacement amplification (dMDA) to generate long molecules with reduced bias. PacBio sequencing generated up to 40% genome coverage per single-cell, enabling detection of single nucleotide variants (SNVs), structural variants (SVs), and tandem repeats, also in regions inaccessible by short reads. 28 somatic SNVs were detected, including one case of mitochondrial heteroplasmy. 5473 high-confidence SVs/cell were discovered, a sixteen-fold increase compared to Illumina-based results from clonally related cells. Single-cell de novo assembly generated a genome size of up to 598 Mb and 1762 (12.8%) complete gene models. In summary, our work shows the promise of long-read sequencing toward characterization of the full spectrum of genetic variation in single cells.

A novel quantitative targeted analysis of X-chromosome inactivation (XCI) using nanopore sequencing.

X-chromosome inactivation (XCI) analyses often assist in diagnostics of X-linked traits, however accurate assessment remains challenging with current methods. We developed a novel strategy using amplification-free Cas9 enrichment and Oxford nanopore technologies sequencing called XCI-ONT, to investigate and rigorously quantify XCI in human androgen receptor gene (AR) and human X-linked retinitis pigmentosa 2 gene (RP2). XCI-ONT measures methylation over 116 CpGs in AR and 58 CpGs in RP2, and separate parental X-chromosomes without PCR bias. We show the usefulness of the XCI-ONT strategy over the PCR-based golden standard XCI technique that only investigates one or two CpGs per gene. The results highlight the limitations of using the golden standard technique when the XCI pattern is partially skewed and the advantages of XCI-ONT to rigorously quantify XCI. This study provides a universal XCI-method on DNA, which is highly valuable in clinical and research framework of X-linked traits.

Transposable element insertions in 1000 Swedish individuals.

The majority of rare diseases are genetic, and regardless of advanced high-throughput genomics-based investigations, 60% of patients remain undiagnosed. A major factor limiting our ability to identify disease-causing alterations is a poor understanding of the morbid and normal human genome. A major genomic contributor of which function and distribution remain largely unstudied are the transposable elements (TE), which constitute 50% of our genome. Here we aim to resolve this knowledge gap and increase the diagnostic yield of rare disease patients investigated with clinical genome sequencing. To this end we characterized TE insertions in 1000 Swedish individuals from the SweGen dataset and 2504 individuals from the 1000 Genomes Project (1KGP), creating seven population-specific TE insertion databases. Of note, 66% of TE insertions in SweGen were present at >1% in the 1KGP databases, proving that most insertions are common across populations. Focusing on the rare TE insertions, we show that even though ~0.7% of those insertions affect protein coding genes, they rarely affect known disease casing genes (<0.1%). Finally, we applied a TE insertion identification workflow on two clinical cases where disease causing TE insertions were suspected and could verify the presence of pathogenic TE insertions in both. Altogether we demonstrate the importance of TE insertion detection and highlight possible clinical implications in rare disease diagnostics.

Gustavson syndrome is caused by an in-frame deletion in RBMX associated with potentially disturbed SH3 domain interactions.

RNA binding motif protein X-linked (RBMX) encodes the heterogeneous nuclear ribonucleoprotein G (hnRNP G) that regulates splicing, sister chromatid cohesion and genome stability. RBMX knock down experiments in various model organisms highlight the gene's importance for brain development. Deletion of the RGG/RG motif in hnRNP G has previously been associated with Shashi syndrome, however involvement of other hnRNP G domains in intellectual disability remain unknown. In the current study, we present the underlying genetic and molecular cause of Gustavson syndrome. Gustavson syndrome was first reported in 1993 in a large Swedish five-generation family presented with profound X-linked intellectual disability and an early death. Extensive genomic analyses of the family revealed hemizygosity for a novel in-frame deletion in RBMX in affected individuals (NM_002139.4; c.484_486del, p.(Pro162del)). Carrier females were asymptomatic and presented with skewed X-chromosome inactivation, indicating silencing of the pathogenic allele. Affected individuals presented minor phenotypic overlap with Shashi syndrome, indicating a different disease-causing mechanism. Investigation of the variant effect in a neuronal cell line (SH-SY5Y) revealed differentially expressed genes enriched for transcription factors involved in RNA polymerase II transcription. Prediction tools and a fluorescence polarization assay imply a novel SH3-binding motif of hnRNP G, and potentially a reduced affinity to SH3 domains caused by the deletion. In conclusion, we present a novel in-frame deletion in RBMX segregating with Gustavson syndrome, leading to disturbed RNA polymerase II transcription, and potentially reduced SH3 binding. The results indicate that disruption of different protein domains affects the severity of RBMX-associated intellectual disabilities.

Novel cancer gene discovery using a forward genetic screen in RCAS-PDGFB-driven gliomas.

BACKGROUND: Malignant gliomas, the most common malignant brain tumors in adults, represent a heterogeneous group of diseases with poor prognosis. Retroviruses can cause permanent genetic alterations that modify genes close to the viral integration site. METHODS: Here we describe the use of a high-throughput pipeline coupled to the commonly used tissue-specific retroviral RCAS-TVA mouse tumor model system. Utilizing next-generation sequencing, we show that retroviral integration sites can be reproducibly detected in malignant stem cell lines generated from RCAS-PDGFB-driven glioma biopsies. RESULTS: A large fraction of common integration sites contained genes that have been dysregulated or misexpressed in glioma. Others overlapped with loci identified in previous glioma-related forward genetic screens, but several novel putative cancer-causing genes were also found. Integrating retroviral tagging and clinical data, Ppfibp1 was highlighted as a frequently tagged novel glioma-causing gene. Retroviral integrations into the locus resulted in Ppfibp1 upregulation, and Ppfibp1-tagged cells generated tumors with shorter latency on orthotopic transplantation. In human gliomas, increased PPFIBP1 expression was significantly linked to poor prognosis and PDGF treatment resistance. CONCLUSIONS: Altogether, the current study has demonstrated a novel approach to tagging glioma genes via forward genetics, validating previous results, and identifying PPFIBP1 as a putative oncogene in gliomagenesis.

Migrating to Long-Read Sequencing for Clinical Routine BCR-ABL1 TKI Resistance Mutation Screening.

Objective: The aim of this project was to implement long-read sequencing for BCR-ABL1 TKI resistance mutation screening in a clinical setting for patients undergoing treatment for chronic myeloid leukemia. Materials and Methods: Processes were established for registering and transferring samples from the clinic to an academic sequencing facility for long-read sequencing. An automated analysis pipeline for detecting mutations was established, and an information system was implemented comprising features for data management, analysis and visualization. Clinical validation was performed by identifying BCR-ABL1 TKI resistance mutations by Sanger and long-read sequencing in parallel. The developed software is available as open source via GitHub at https://github.com/pharmbio/clamp. Results: The information system enabled traceable transfer of samples from the clinic to the sequencing facility, robust and automated analysis of the long-read sequence data, and communication of results from sequence analysis in a reporting format that could be easily interpreted and acted upon by clinical experts. In a validation study, all 17 resistance mutations found by Sanger sequencing were also detected by long-read sequencing. An additional 16 mutations were found only by long-read sequencing, all of them with frequencies below the limit of detection for Sanger sequencing. The clonal distributions of co-existing mutations were automatically resolved through the long-read data analysis. After the implementation and validation, the clinical laboratory switched their routine protocol from using Sanger to long-read sequencing for this application. Conclusions: Long-read sequencing delivers results with higher sensitivity compared to Sanger sequencing and enables earlier detection of emerging TKI resistance mutations. The developed processes, analysis workflow, and software components lower barriers for adoption and could be extended to other applications.

Loss of Nexilin function leads to a recessive lethal fetal cardiomyopathy characterized by cardiomegaly and endocardial fibroelastosis.

The Nexilin F-Actin Binding Protein (Nexilin) encoded by NEXN is a cardiac Z-disc protein important for cardiac function and development in humans, zebrafish, and mice. Heterozygote variants in the human NEXN gene have been reported to cause dilated and hypertrophic cardiomyopathy. Homozygous variants in NEXN cause a lethal form of human fetal cardiomyopathy, only described in two patients before. In a Swedish, four-generation, non-consanguineous family comprising 42 individuals, one female had three consecutive pregnancies with intrauterine fetal deaths caused by a lethal form of dilated cardiomyopathy. Whole-exome sequencing and variant analysis revealed that the affected fetuses were homozygous for a NEXN variant (NM_144573:c.1302del;p.(Ile435Serfs*3)). Moreover, autopsy and histology staining declared that they presented with cardiomegaly and endocardial fibroelastosis. Immunohistochemistry staining for Nexilin in the affected fetuses revealed reduced antibody staining and loss of striation in the heart, supporting loss of Nexilin function. Clinical examination of seven heterozygote carriers confirmed dilated cardiomyopathy (two individuals), other cardiac findings (three individuals), or no cardiac deviations (two individuals), indicating incomplete penetrance or age-dependent expression of dilated cardiomyopathy. RNA sequencing spanning the variant in cDNA blood of heterozygote individuals revealed nonsense-mediated mRNA decay of the mutated transcripts. In the current study, we present the first natural course of the recessively inherited lethal form of human fetal cardiomyopathy caused by loss of Nexilin function. The affected family had uneventful pregnancies until week 23-24, followed by fetal death at week 24-30, characterized by cardiomegaly and endocardial fibroelastosis.

The Value of Whole-Genome Sequencing for Mitochondrial DNA Population Studies: Strategies and Criteria for Extracting High-Quality Mitogenome Haplotypes.

Whole-genome sequencing (WGS) data present a readily available resource for mitochondrial genome (mitogenome) haplotypes that can be utilized for genetics research including population studies. However, the reconstruction of the mitogenome is complicated by nuclear mitochondrial DNA (mtDNA) segments (NUMTs) that co-align with the mtDNA sequences and mimic authentic heteroplasmy. Two minimum variant detection thresholds, 5% and 10%, were assessed for the ability to produce authentic mitogenome haplotypes from a previously generated WGS dataset. Variants associated with NUMTs were detected in the mtDNA alignments for 91 of 917 (~8%) Swedish samples when the 5% frequency threshold was applied. The 413 observed NUMT variants were predominantly detected in two regions (nps 12,612-13,105 and 16,390-16,527), which were consistent with previously documented NUMTs. The number of NUMT variants was reduced by ~97% (400) using a 10% frequency threshold. Furthermore, the 5% frequency data were inconsistent with a platinum-quality mitogenome dataset with respect to observed heteroplasmy. These analyses illustrate that a 10% variant detection threshold may be necessary to ensure the generation of reliable mitogenome haplotypes from WGS data resources.

CRISPR-Cas9 induces large structural variants at on-target and off-target sites in vivo that segregate across generations.

CRISPR-Cas9 genome editing has potential to cure diseases without current treatments, but therapies must be safe. Here we show that CRISPR-Cas9 editing can introduce unintended mutations in vivo, which are passed on to the next generation. By editing fertilized zebrafish eggs using four guide RNAs selected for off-target activity in vitro, followed by long-read sequencing of DNA from >1100 larvae, juvenile and adult fish across two generations, we find that structural variants (SVs), i.e., insertions and deletions ≥50 bp, represent 6% of editing outcomes in founder larvae. These SVs occur both at on-target and off-target sites. Our results also illustrate that adult founder zebrafish are mosaic in their germ cells, and that 26% of their offspring carries an off-target mutation and 9% an SV. Hence, pre-testing for off-target activity and SVs using patient material is advisable in clinical applications, to reduce the risk of unanticipated effects with potentially large implications.

A polygenic risk score predicts mosaic loss of chromosome Y in circulating blood cells.

BACKGROUND: Mosaic loss of Y chromosome (LOY) is the most common somatic change that occurs in circulating white blood cells of older men. LOY in leukocytes is associated with increased risk for all-cause mortality and a range of common disease such as hematological and non-hematological cancer, Alzheimer's disease, and cardiovascular events. Recent genome-wide association studies identified up to 156 germline variants associated with risk of LOY. The objective of this study was to use these variants to calculate a novel polygenic risk score (PRS) for LOY, and to assess the predictive performance of this score in a large independent population of older men. RESULTS: We calculated a PRS for LOY in 5131 men aged 70 years and older. Levels of LOY were estimated using microarrays and validated by whole genome sequencing. After adjusting for covariates, the PRS was a significant predictor of LOY (odds ratio [OR] = 1.74 per standard deviation of the PRS, 95% confidence intervals [CI] 1.62-1.86, p < 0.001). Men in the highest quintile of the PRS distribution had > fivefold higher risk of LOY than the lowest (OR = 5.05, 95% CI 4.05-6.32, p < 0.001). Adding the PRS to a LOY prediction model comprised of age, smoking and alcohol consumption significantly improved prediction (AUC = 0.628 [CI 0.61-0.64] to 0.695 [CI 0.67-0.71], p < 0.001). CONCLUSIONS: Our results suggest that a PRS for LOY could become a useful tool for risk prediction and targeted intervention for common disease in men.

Phylogenetic history of patrilineages rare in northern and eastern Europe from large-scale re-sequencing of human Y-chromosomes.

The most frequent Y-chromosomal (chrY) haplogroups in northern and eastern Europe (NEE) are well-known and thoroughly characterised. Yet a considerable number of men in every population carry rare paternal lineages with estimated frequencies around 5%. So far, limited sample-sizes and insufficient resolution of genotyping have obstructed a truly comprehensive look into the variety of rare paternal lineages segregating within populations and potential signals of population history that such lineages might convey. Here we harness the power of massive re-sequencing of human Y chromosomes to identify previously unknown population-specific clusters among rare paternal lineages in NEE. We construct dated phylogenies for haplogroups E2-M215, J2-M172, G-M201 and Q-M242 on the basis of 421 (of them 282 novel) high-coverage chrY sequences collected from large-scale databases focusing on populations of NEE. Within these otherwise rare haplogroups we disclose lineages that began to radiate ~1-3 thousand years ago in Estonia and Sweden and reveal male phylogenetic patterns testifying of comparatively recent local demographic expansions. Conversely, haplogroup Q lineages bear evidence of ancient Siberian influence lingering in the modern paternal gene pool of northern Europe. We assess the possible direction of influx of ancestral carriers for some of these male lineages. In addition, we demonstrate the congruency of paternal haplogroup composition of our dataset with two independent population-based cohorts from Estonia and Sweden.

Immune cells lacking Y chromosome show dysregulation of autosomal gene expression.

Epidemiological investigations show that mosaic loss of chromosome Y (LOY) in leukocytes is associated with earlier mortality and morbidity from many diseases in men. LOY is the most common acquired mutation and is associated with aberrant clonal expansion of cells, yet it remains unclear whether this mosaicism exerts a direct physiological effect. We studied DNA and RNA from leukocytes in sorted- and single-cells in vivo and in vitro. DNA analyses of sorted cells showed that men diagnosed with Alzheimer's disease was primarily affected with LOY in NK cells whereas prostate cancer patients more frequently displayed LOY in CD4 + T cells and granulocytes. Moreover, bulk and single-cell RNA sequencing in leukocytes allowed scoring of LOY from mRNA data and confirmed considerable variation in the rate of LOY across individuals and cell types. LOY-associated transcriptional effect (LATE) was observed in ~ 500 autosomal genes showing dysregulation in leukocytes with LOY. The fraction of LATE genes within specific cell types was substantially larger than the fraction of LATE genes shared between different subsets of leukocytes, suggesting that LOY might have pleiotropic effects. LATE genes are involved in immune functions but also encode proteins with roles in other diverse biological processes. Our findings highlight a surprisingly broad role for chromosome Y, challenging the view of it as a "genetic wasteland", and support the hypothesis that altered immune function in leukocytes could be a mechanism linking LOY to increased risk for disease.

Characterization of the nuclear and cytosolic transcriptomes in human brain tissue reveals new insights into the subcellular distribution of RNA transcripts.

Transcriptome analysis has mainly relied on analyzing RNA sequencing data from whole cells, overlooking the impact of subcellular RNA localization and its influence on our understanding of gene function, and interpretation of gene expression signatures in cells. Here, we separated cytosolic and nuclear RNA from human fetal and adult brain samples and performed a comprehensive analysis of cytosolic and nuclear transcriptomes. There are significant differences in RNA expression for protein-coding and lncRNA genes between cytosol and nucleus. We show that transcripts encoding the nuclear-encoded mitochondrial proteins are significantly enriched in the cytosol compared to the rest of protein-coding genes. Differential expression analysis between fetal and adult frontal cortex show that results obtained from the cytosolic RNA differ from results using nuclear RNA both at the level of transcript types and the number of differentially expressed genes. Our data provide a resource for the subcellular localization of thousands of RNA transcripts in the human brain and highlight differences in using the cytosolic or the nuclear transcriptomes for expression analysis.